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SHIELDOSE: A Computer Code for Space-Shielding Radiation Dose Calculations! 


Stephen M. Seltzer 
National Bureau of Standards 
Washington, D.C. 20234 


A computer code, SHIELDOSE, has been developed for the 
calculation of absorbed dose as a function of depth in aluminum 
shielding material of spacecraft, given the electron and proton 
fluences encountered in orbit. Absorbed dose, for small volumes 
of the detector materials Ag, HA0, Si, and Si0,, is evaluated in 
three geometries: (1) in a semi-infinite plane medium, (2) at 
the transmission surface of a finite-thickness slab, and (3) at 
the center of a solid sphere. Use of the code is described, and 
an extensive set of monoenergetic depth-dose data for the various 
detector materials and geometries is tabulated. 


Key words: Computer code; depth-dose data; electrons; electron 
bremsstrahlung; protons; space shielding. 


1. Introduction 
The ability to predict absorbed dose within a spacecraft due to a 
specified radiation environment is important for design and planning consid- 
erations pertaining to the reliability of electronic components and to the 
radiological safety of on-board personnel. This report describes the computer 
code SHIELDOSE, which evaluates the absorbed dose as a function of depth in 
aluminum shielding material of spacecraft, given the electron and proton 


fluences encountered in orbit. 


Mork Supported by the Office of Naval Research and by the National Space 
Science Data Center at the NASA Goddard Space Flight Center. 
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The code makes use of pre-calculated, monoenergetic depth-dose data for 
an isotropic, broad-beam fluence of radiation incident on uniform aluminum 
plane media. Such data are particularly suitable for routine dose predictions 
in situations where the geometrical and compositional complexities of the 
spacecraft are not known. Furthermore, restricting our consideration to these 
rather simple geometries has allowed for the development of accurate electron 
and electron-bremsstrahlung data sets based on detailed transport calculations 
rather than on more approximate methods. The present version of SHIELDOSE 
calculates, for arbitrary proton and electron incident spectra, the dose 
absorbed in small volumes of the detector materials Ag, Ho0 (tissue-equivalent 
detector), Si, and Si0,, in the following aluminum shield geometries: (1) in 
a semi-infinite plane medium, as a function of depth; (2) at the transmission 
surface of a plane slab, as a function of slab thickness; and (3) at the 
center of a solid sphere, as a function of sphere radius. 

Section 2 defines a number of pertinent quantities used in the calcula- 
tions and briefly outlines the methods used for integrating the monoenergetic 
data over the incident flux spectra. Section 3 gives the meanings and 
required formats of the various input parameters, and other running 
instructions for SHIELDOSE. Appendix A is a reprint of a report fine which 
gives details on the methods with which the monoenergetic depth-dose 
distributions were calculated. Appendix B comprises tables of the mono- 
energetic depth-dose data; appendix C gives the input and output listings of a 


sample run; and appendix D contains a FORTRAN listing of the code. 


2Figures in brackets indicate the literature references at the end of this 
paper. 


2. Method 
The quantity of interest calculated by SHIELDOSE is the absorbed dose 
9(z), as a function of depth z in an aluminum shield, for an arbitrary 


incident fluence? of radiation. The incident fluence can be either protons or 


' electrons. If it is protons, it will be implicitly understood that 9(z) is 


the proton dose; if electrons, 9(z) is calculated separately for the 


bremsstrahlung component and the "electron" component (by which we mean the 


| distribution of deposited electron energy never converted to bremsstrahlung 


Case l. Q 


photons). Results are calculated for three geometries (see fig. 1): 


det (7) is the absorbed dose in a thin detector (of material 


co 


denoted by "det") at a depth z ina semi-infinite 
plane aluminium medium. Irradiation is from one side 
only (the assumed infinite backing effectively insures 


this). 


Case 2. gt (2) is the absorbed dose in a thin detector at the trans- 
mission surface of an aluminum plane slab of thickness 


z. Irradiation is from one side only. 


Case 3. gt (7) is the absorbed dose in a small detector volume at the 
center of a solid aluminum sphere of radius Zz. 


Irradiation is from all directions. 


The detector materials AX, H,0, Si, and Si0, are presently. included in the 


data base. 


3We can use the fluence (number per unit area, per unit energy, per unit 
solid angle) and the flux (number per unit area, per unit energy, per unit 
solid angle, per unit time) interchangeably. Use of the fluence.leads to 
absorbed dose, and the flux to dose rate. 


We define the following quantities: 


o(T, 28) is the incident fluence, differential in energy We and 


2 -] 


angle 6, having units of cm Mev! sr We assume 


to be isotropic and drop the explicit dependence on 6. 


aq iT ) = 40 o(T,) is the full-space (47) incident fluence, having units 


re) 
of cm~2Mev7! 


Doe za) is the depth-dose distribution for irradiation by a 


foe) 


monoenergetic, omnidirectional beam with energy Ty 
incident on a semi-infinite medium (case 1 geometry). 
The depth-dose distribution, normalized to an incident 
current of 1 m7? (see appendix A, eqs (1-5)), has units 


of MeV om@/g. Daa 


co 


is given separately for the proton, 


electron, and electron-bremsstrahlung components. 


po (2.0) is similar to pore but for monoenergetic irradiation 
of finite-thickness slabs (case 2 geometry), from one 
side only. 
Then 
det 2K det 
afet(2) = Ff vaglty) O8F2.7,) aT, (1) 


where K = 1.6022 x 107° rads/ (MeV any and the factor of 1/4 comes from the 
relation between the full-space fluence and the incident current (see eq (5) 
in appendix A). Because the incident fluence Vag (Tg) typically falls off 


rather steeply with increasing energy, eq (1) is rewritten 


Qn T 
i aa Teepe yebeee ei, dltngligh (2) 
An min ei 


P|~ 


| Equation (2) is evaluated numerically using a Simpson's rule, so that 


N 
gi*t(z) = £4 De Wig Ygq(Ty) DERN ZT:) (3) 
i=1. 


where w, = 1,4,2,4,2 ... 2,41 (Nodd),” and gn T, = 2nT + (i - 1)A 
/(N- 1). 


In order to ee the interpolation of D 


with A = (n i mi) 


ger 


co 


Z5T,) among 


tabulated values, eq (3) is recast as 


N 
giz) - £4 2 miTy YanlTy) 8 SNT,) {8(T,) DE *Ex(z.7)). Typ » (4) 


for which the scaling functions are 


S(T) : x(z,T)_ 
protons A/a eg cm7? Mev7 z/¥v(T) 
electrons PM af 9 cm™* Mev! 2/5) 
bremsstrahlung /tigexs Mey"! z/T 


and where i is the mean range evaluated in the continuous-slowing-down 
approximation. Values of the mean range Gay ane given in table 1 for 
protons [2] and table 2 for electrons [3]. The scaled depth-dose distri- 


butions in the curly brackets in eq (4), having a much reduced dependence 


4 any N22 can be handled, and the weights w. 


; are suitably adjusted. 


On z and Ty: are tabulated in appendix B. Interpolation is accomplished 
by fitting with natural cubic splines [4], first in T, then in x. 


Yalues for T and Teas (eq (2) aboye) are often imposed by the 


min 
available knowledge of the incident electron or proton fluence Van (To)> cfs 

in the evaluation of eq (4), T; is outside of the energy range for which the 
scaled, monoenergetic depth-dose distributions are giyen, then the curly 

bracket is evaluated for the nearest energy tabulated, and the energy 

variation is assumed to be taken care of by the term Sin): This procedure 
should cause little error since the tabulated proton data extend from 2 to 

5000 MeV, the electron data from 0.1 to 10 MeV, and the bremsstrahlung data 

from 0.02 to 20 MeV, and should cover incident spectra encountered in space- 
shielding problems. 

A more likely situation is that information on the incident flux spectra 
extends only over an energy region more limited than those above. For example, 
current model predictions of the terrestrial radiation environment [5] typically 
give spectral information over the energy ranges 0.] to 10 MeV for electrons 
0.5 to 500 MeV for trapped protons, and 10 to 200 MeV for solar-flare protons. 
If qq (To) is entered into SHIELDOSE in numerical form, care should be 
taken in choosing values of fn and Tray in order to avoid the automatic 
cubic-spline extrapolation of Van (To) which can occasionally result in 
Significant, and quite peculiar, components in the incident energy spectrum. 

It should be kept in mind that the final results neglect those contributions 


to the dose from radiation with Us <al! and Ue > ax This can be 


min 


particularly important for depths z<r_(T_:._). Some results of test 


Oo’ min 
calculations using 


are given in table 3 which shows, at various depths z and for various values 


of the e-folding energy a, the effect on g*(z) of different choices of 
Tin and Tax’ 
The discussion above outlines the calculation of ecu) for a semi- 


foe) 


infinite medium. The calculation of @ for finite-thickness slabs is 


pdet 


det (7) 


similar but based on (2,15). The results for the dose at the center of a 


solid sphere of radius z are obtained from (see also eq (10) in appendix A) 


d log g*t(2') 
muau(aancgec’(z) 4) = ————— (5) 
eo d log z' 


z'=2 


where the logarthmic derivative is evaluated using a cubic-spline fit of 


log gt (2), We have found that N = 301 for protons and N = 101 for 
electrons gives distributions of gt(7) sufficiently smooth to give 
reasonably smooth results for git(2). 


3. Input Parameters 


The format of the input deck is (using FORTRAN notation): 


IDET, NUT, IPRNT 316 


[scaled monoenergetic depth-dose data set]” not user supplied 


IMAX, IUNT 216 
(Z(1), I=1, IMAX) 6F12.5 
EMINS, EMAXS, EMINP, EMAXP, NPTSP, AFNORS Os meralOn sro 


EMINE, EMAXE, NPTSE 


“This rather lengthy data set, about 13,400 words of information, should be 


permanently stored on tape or disk. Although its position in the input 
Stream is as indicated above, it is read from unit NUT. 


The following block of input data can be repeated for each set of incident 


spectra: 

COMMENT 72H 

JSMAX, JPMAX, JEMAX, EUNIT, TINTER S16, -IPZENZe0 
(EPSS(J), J=1, JSMAX) 1P6E12.4 
(SS(J), J=1, JSMAX) 1P6E12.4 
(EPSP(J), J=1, JPMAX) 1P6E12.4 
(SP(J), J=1, JPMAX) 1P6E12.4 
(ERSE()5 » JEMAX) 1P6E12.4 
(SE(J), J=1, JEMAX) 1P6E12.4 


The meaning of the input parameters is as follows: 
IDET = 1 , detector material is AX. 
2 , detector material is H50. 
3, detector material is Si. 


4 , detector material is Si0,. 


NUT is the unit from which the scaled depth-dose data set is 


read. 


IPRNT = 1 , input monoenergetic depth-dose data is listed in 
printout. 
2 , input monoenergetic depth-dose data is not listed in 


printout. 


IMAX is the number of depths z _ for which dose calculation is 


desired (IMAX < 50). 


IUNT = 


EMINS 


EMAXS 


EMINP 


EMAXP 


NPTSP 


lease GeDthe Zain mils.° 
ar, eaeptn 2 >in g/cm*. 


3, depth z inm. 


is simultaneously the IMAX depths (case 1 geometry), 
thicknesses (case 2 geometry), and sphere radii (case 3 
geometry) in units specified by IUNT. No Z-value 

can be O because of a log Z_ interpolation performed 
in the bremsstrahlung calculation. Furthermore, 
0<Z<«< 30 g/cm*, due to the extent of the bremsstrah- 
lung data. 


is the lower limit T for the integration over the 


min 
incident solar-flare-proton spectrum (see eq (2)). 


is the upper limit ran for the integration over the 


incident solar-flare-proton spectrum (see eq (2)). 


is the lower limit ae for the integration over the 


incident trapped-proton spectrum (see eq (2)). 


is the upper limit Uae for the integration over the 


incident trapped-proton spectrum (see eq (2)). 


is the number of points used in the numerical evaluation 


of the integral over the incident proton spectrum, for 


0.001 inch =,0.0254 mm. For aluminum, with a density of 2.70 g/cm? 


0.00686 g/cm. 


2 


NPTSE 


JPMAX 


which the integration extends from min{EMINS. DMINP) to 
max(EMAXS, EMAXP). Note that EMINS, EMAXS. EMINP, and 
EMAXP may be slightly adjusted to the nearest mesh- 
energy values. (NPTSP < 30], recommend WPTSP = 301). 
for the integration over the 


as the jower limit 1 ie 


incident electron spectrum (see eg (2))- 


is the upper limit i‘. for the integration over the 
incident electron spectrum (see eg (2))-. 


is the number of points wsed in the numerical evaluation 
of the integral over the incident electron spectrum 
(NPTSE < 7203, recommend WNPTSE = 10])- 


is any identifying message. 


is the number of energies for which the incident solar— 
flare-proton spectrum is read in (JSMAK< 10). If 
JSMAX <2, mo sOlar-flare-proton results are celculatec. 


is the number of energies for which the incident trapped- 
proton spectrum is read in (JPMAK< 101). If 
JPMAX < Z, mo trapped-proton results are calculated. 


is the number of energies for which the incident elecdron 
spectrum is read in (JEMAK< i101). If JEMAK <2. mm 
electron and bremsstrahlung results are calculated. 
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EUNIT 


TINTER 


EPSS 


is the conversion factor from (unit energy)” to Mey"! , 


where (unit energy)! pertains to the units in which 
the incident proton and electron flux spectra are read. 
For example, if the incident spectra are read in units of 


key! | 


, then EUNIT = 1000; if read in units of MeV 
then EUNIT = 1 (if the EUNIT field is left blank or if 


EUNIT <0, then SHIELDOSE assumes EUNIT = 1). 


is the time interval (e.g., mission duration) for which 
results are desired. TINTER is in units of (unit time), 
where (unit time) 7! pertains to the units in which the 
incident trapped-proton and electron flux spectra are 
read. For example, if incident spectra are read in units 
of sec”! , then TINTER = 1 gives the absorbed dose for 
a l-sec irradiation; TINTER = 3.1536 x 10’ gives dose 
for one year (365 days); etc. If the TINTER field is 
left blank or if TINTER <0, then SHIELDOSE assumes 
TINTER = 1. Note that TINTER does not pertain to the 
incident solar-flare-proton fluence spectrum which is 
assumed to be the result of a statistical analysis 
relating the probability of the incident solar-flare- 


proton fluence to the mission duration [6]. 
are the JSMAX energies in MeV (in ascending order) for 


which the incident solar-flare-proton fluence spectrum is 


read. If JSMAX < 2, this array is not read in. 


11 


SS 


EPSP 


SP 


EPSE 


SE 


are the values (at energies EPSS) of the full-space 
(41), omnidirectional solar-flare-proton incident fluence 
spectrum in units of cm72 (unit energy) 7). 

are the JPMAX energies in MeV (in ascending order) for 
which the incident trapped-proton flux spectrum is read. 


If JPMAX < 2, this array is not read in. 


are the values (at energies EPSP) of the full-space 
(47), omnidirectional trapped-proton incident flux 


spectrum in units of cm72 (unit energy)” (unit time)~!. 


are the JEMAX energies in MeV (in ascending order) 


for which the incident electron fluence spectrum is read. 


are the values (at energies EPSE) of the full-space 
(47), omnidirectional incident electron fluence spectrum 


in units of m7? (unit energy) 7! (unit time)~!. 


Note the following option which pertains to the solar- 
flare-proton, trapped-proton, and electron incident 


spectra: 


If EPS(1) <0 (but JMAX > 3), then the incident 
spectrum is assumed to be Van (To) = (C/a) exp(- T,/o)» 
with o(MeV) = S(1) and C = S(2). The units for C 


are cm7* for the solar-flare-proton fluence spectrum, 


and cm” (unit time)! for the trapped-proton and 
electron flux spectra. This option may be particularly 
useful for the solar-flare-proton spectrum for which, in 
the active years 1977-1983, the spectrum is characterized 
by a= 26.5 MeV for anomalously large events and by 

a = 63.3 MeV for ordinary events [6]. Values of C can 
be obtained from information on the corresponding integral 
spectrum I(> To): via C = I(> i) exp(T,/a). If the 
entered value S(2) <0, then SHIELDOSE sets C= 1. 


The code is written to process quickly many incident spectra in a single 
run. For example, using an IBM 360/91 computer and with NPTSP = 301, 

NPTSE = 101 and IMAX = 48, approximately 10 sec is required for setting up 
internal matrices and ~ 1] sec for each input proton and electron spectrum 
set. The storage requirement is ~ 80 K words. 

It is, in general, necessary to suppress underflow diagnostics which 
interrupt, and can eventually terminate, a run. Underflow conditions are not 
detected by the code itself and should therefore be suppressed using the 
appropriate options available at the particular computing installation. For 
the IBM system used by the author, this is accomplished by the "CALL ERRSET" 


instruction at the beginning of the main routine listed in appendix D. 


REFERENCES 
[1] Seltzer, S. M., IEEE Trans. Nucl. Science, NS-26, 4896 (1979). 
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79-01; all from the National Space Science Data Center, Greenbelt, 
Maryland. 
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Table 7. Proton range mu) as 2 function of protom energy T,- 
Values are from reference [2]. 


T, (Mev) (afar) Tv) sr, (a/em”) 
2 1.147 x 10% 100 9.975 
3 2.181 x 107 150 2.026 x 10° 
4 3.500 x 107 200 3.320 x to! 
5 5.094 x 10-7 300 6.547 x To! 
6 6.574 x 10 - 400 1.043 x 10” 
g 1.157 x 10! 500 1.482 x 10° 
10 1.703 x 10°! =) 1.959 x 107 
15 3.458 x 10 | g00 2.994 x 10" 
20 5.739 x 10! 1000 4.102 x 10" 
30 1.177 1500 7.046 x 10- 
7) 1.963 2000 1.010 x 16° 
= 2.920 3000 1.627 x 10° 
60 4.039 4000 2.239 x 10° 
50 6.728 5006 2.943 x 10° 


Table 2. Electron range ly as a function of electron energy ee 


Values are from reference [3]. 


T, (MeV) r(g/em*) T, (MeV) r(g/em®) 
0.01 3.519 x 107" 1.0 5.493 x 107 
0.015 7.074 x 1074 1.5 8.825 x 107! 
3.02 1.165 x 107° 2.0 Teale 
0.03 2.356 x 107° 3.0 1.855 
0.04 3.883 x 107° 4.0 2.476 
0.05 5.714 x 107° Qt Sue 
0.06 7822) AO 6.0 3.658 
0.08 1,279 x 107° 8.0 4.777 
0.10 1.864 x 1072 10 5.84] 

0.15 3.641 x 107° 15 8.305 
0.20 5772 x 1084 20 1.054 x 10! 
0.30 1.077 x 107! 30 1.448 x 10! 
0.40 1.640 x 107! 40 1.788 x 10! 
0.50 2.243 x 107! 50 2.087 x 10! 
0.60 2.871 x 107! 60 2.355 x 10! 
0.80 4.168 x 107! 80 2.816 x 10! 
100 3.204 x 10! 
16 


Fraction of dose a (2) contributed by incident spectral energies below 


1: The quantity given corresponds to the ratio 


T [o0] 

c pAL x i Ag 

i ap (z,T,) exp ( T,/o) dT / I De (z,T)) exp(-T,/a) qT ; 

and is given for various depths z in g/cm? and for various a and 
T. in MeV. For any z and a, the fraction of the dose contributed 
Dye = T= 7 is the difference between the value tabulated for 
= 7, and that for Wes =Paly 


(e max 
electrons, and bremsstrahlung. 


min: Results are given for protons, 


PROTONS 
2 5 10 200 500 1000 2000 5000 
a = 10 MeV 
0.00 0.15 0.63 1.00 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 1.00 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 1.00 100 1.00 1.00 1.00 
0.00 0.00 0.00 0.51 0.92 0.96 0.99 1.00 
a = 20 MeV 
0.00 0.10 0.47 1.00 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 1.00 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 1.00 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 0.27 0.97 0.95 0.99 1.00 
a = 50 MeV 
0.00 0.06 0.31 1.00 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 0.99 1.00 1.00 100 1.00 
0.00 0.00 0.00 0.97 1.00 1.00 1600 1.00 
0.00 0.00 0.00 0.10 0.99 1.00 1.00 1.00 
o = 100 MeV 
0.00 0.05 0.24 0.96 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 0.95 1.00 1.00 1.00 1.00 
0.00 0.00 0.00 0.86 0.99 1.00 1.00 1.00 
0.00 0.00 0.00 0.05 0.95 1.00 1.00 1.00 
a = 200 MeV 
0.00 70503 0.18 0.87 0.98 1.00 1.00 1.00 
0.00 0.00 0.00 0.83 0.97 1.00 1.00 1.00 
0.00 0.00 0.00 0.69 0.95 1.00 1.00 1.00 
0.00 0.00 0.00 0.03 0.81 0.99 1.00 1.00 
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Table 3. (Continued) 

ELECTRONS 
0.5 2 

a = 0.25 MeV 
OA 1.00 
0.00 0.98 
0.00 0.00 

a = 0.50 MeV 
0.48 0.98 
0.00 0.84 
0.00 0.00 

a = 1 Mey 
0.29 0.85 
0.00 O57: 
0.00 0.00 

a = 2 MeV 
0.16 0.62 
0.00 0.33 
0.00 0.00 

a = 4 MeV 
0.09 0.39 
0.00 (OI 7/ 
0.00 0.00 
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Table 3. (Continued) 
BREMSSTRAHLUNG 
0.5 2 
= 0.25 MeV 
0.68 1.00 
O33 0.99 
0.30 0.98 
0.16 0.97 
= 0.50 MeV 
0.40 0.96 
0.09 0.85 
{0)-<(0)7/ 0.74 
0.03 0.62 
= ] Mey 
0.19 0.74 
0.02 0.51 
0.01 0.29 
0.00 0.18 
= 2 MeV 
0.08 0.40 
0.01 0.22 
0.00 0.07 
0.00 0.03 
= 4 MeV 
0.02 OF 15 
0.00 0.08 
0.00 0.01 
0.00 0.00 
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INCIDENT PROTONS 
AND ELECTRONS 


INCIDENT PROTONS 
AND ELECTRONS 


DETECTOR 


CaseE 2. FINITE-THICKNESS SLAB; ABSORBED DOSE IN DETECTOR DENOTED AS 2°t (7) 


INCIDENT PROTONS 
AND ELECTRONS 


det 
Case 3, SOLID SPHERE; ABSORBED DOSE IN DETECTOR DENOTED AS Do (z) 


Fig. 1. Geometries considered in SHIELDOSE calculations. Isotropic, 
broad-beam fluxes of protons and electrons are assumed incident 
on aluminum targets; absorbed dose is calculated for small 
volumes of detector materials AX, H,0, Si, and Si0,. 
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ELECTRON, ELECTROWN-BREMSSTRAHLUNG AND PROTON DEPTH-DOSE 
DATA FOR SPACE-SHIELDING APPLICATIONS* 


Stephen M. Seltzer 
National Bureau of Standards 
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Abstract 


A data set has been developed, consisting of 
depth-dose distributions for omni-directional electron 
and proton fluxes incident on aluminum shields. The 
principal new feature of this work is the accurate 
treatment, based on detailed Monte Carlo calculations, 
of the electron-produced bremsstrahlung component. 
Results covering the energy region of interest in 
space-shielding calculations have been obtained for the 
absorbed dose (a) as a function of depth in a semi- 
infinite medium, (b) at the edge of slab shields, and 
(c) at the center of a solid sphere. The dose to a 
thin tissue-equivalent detector was obtained as well as 
that in aluminum. Various results and comparisons with 
other work are given. 


Introduction 


The longer durations of currently planned and 
future satellite missions, and the increasing preference 
for the use of radiation-sensitive CMOS (Complimentary 
Metal-Oxide-Semiconductor) devices, make more acute the 
problem of predicting the radiation dose at points in 
the spacecraft.! Such predictions. depend on 
(a) knowledge of the spectra of electrons and protons 
incident on the spacecraft; (b) details of the geomet- 
rical and compositional complexities of the spacecraft; 
and (c) an accurate description of the transport of the 
radiation through the spacecraft material. 


Improvements in space-shielding calculations are 
desirable in all three of these areas. The uncertainty 
in the information of the Earth's radiation environment 
is presently approaching a level of about a factor of 
two, and should improve further as more comprehensive 
experimental data are accumulated and theoretical 
models are refined. The complexities of the spacecraft 
configuration can, in principle, be taken into account 
by Monte Carlo transport calculations in complicated 
geometry* or by irradiations of satellite mock-ups. 
However, a far less costly alternative is the use of 
approximate methods*® based on data obtained for simple 
geometries. 


In this paper, calculated results are presented on 
the penetration of electrons and protons in aluminum, 
as summarized in terms of depth-dose data, for condi- 
tions pertinent to space-shielding problems. The 
principal new feature of this work is the accurate 
treatment of the production and transport of electron- 
produced bremsstrahlung photons which are more penetra- 
ting than the electrons themselves. Information on 
the bremsstrahlung tail can be particularly important 
for geostationary orbits in which satellites encounter 
few or no trapped protons, and for which the dose 
penetrating through thick shields is completely due to 
the bremsstrahlung. 


* 

Work supported by the Office of Naval Research and the 
National Space Science Data Center at the NASA Goddard 
Space Flight Center. 


Although the basic calculations were done for 
semi-infinite aluminum targets, consideration has been 
given also to other simple configurations, namely the 
dose at the edge of finite aluminum slabs and the dose 
at the center of uniform, solid aluminum spheres. In 
addition, the dose has been calculated in thin tissue- 
equivalent detectors embedded in or at the edge of the 
aluminum target. 


Method of Calculation 


Electrons and Bremsstrahlung 


The electron calculations make use of the Monte 
Carlo code ETRAN which treats the following processes: 
i) electron energy loss, including energy-loss straggling 
(fluctuations) due both to multiple inelastic scattering 
by atomic electrons and to the emission of bremsstrahlung 
photons; ii) angular deflections of electrons due to 
multiple elastic scattering by atoms; iii) penetration 
and diffusion of the secondary bremsstrahlung photons; 
and iv) penetration and diffusion of energetic secondary 
electrons produced in electron-electron knock-on 
collisions (delta rays) and in the interaction of 
bremsstrahlung photons with the medium (pair, Compton, 
and photoelectrons). 


Details of the Monte Carlo model and of the cross 
sections used have been given elsewhere’~’, along with 
numerous comparisons with experimental results. These 
will not be discussed here, except for the following 
improvement. The ETRAN bremsstrahlung results, based 
on the use of a set of empirically-corrected Bethe- 
Heitler bremsstrahlung cross sections®, were adjusted 
to reflect the exact calculations of the bremsstrahlung 
production cross section of Pratt, Tseng, et al®. 
Comparisons, using both the ETRAN and Pratt cross- 
section data sets for aluminum, of the electron mean 
energy loss per unit pathlength due to the emission of 
bremsstrahlung (see, e.g., reference 11) and of the 
total energy radiated in the course of slowing down (in 
the continuous-slowing-down approximation) give a 
multiplicative correction factor for the ETRAN results 
that is very close to unity for electron source energies 
Ty = 1 MeV, and is 0.96, 0.94, 0.92, 0.89, and 0.82 at 
0.5, 0.2, 0.1, 0.05, and 0.02 MeV, respectively. 


Protons 


The main emphasis on this work was the calculation 
of the dose from electrons. However, it was considered 
desirable to include some corresponding results for 
protons so that one can put the contributions from the 
two types of radiation into proper perspective. For 
this purpose, we limited ourselves to a simple treatment 
that included Coulomb interactions but neglected nuclear 
interactions. The error incurred by this simplification 
is generally no greater than 10-20% for shields up to 
~ 30 g/cm?2, 


The proton calculations were done in the straight- 
ahead, continuous-slowing-down approximation using the 
stopping-power and range data of Barkas and Berger. ° 
Alsmiller et al!° have shown that the neglect of angular 
deflections and range straggling is negligible in 
space-shielding calculations. 
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Initial and Boundary Conditions 


The calculations were made for radiation whose 
fluence is isotropic and which is incident on a semi- 
infinite, uniform aluminum plane medium. The radiation 
is assumed to be incident over a large area, so that the 
transport problem reduces to a one-dimensional one, with 
the depth z as the only spatial variable. Radiation 
crossing out through the plane boundary (z= 0) was 
considered lost to the medium (backscattered). 


Let (Ty,8) be the incident fluence, differential 
in initial energy Tg and in angle €. The fiuence, 
having units of cm “MeV ‘sr *, pertains to the number 
of particles crossing through a unit area perpendicular 
to the incident particle vector. The current of 
particles crossing into the medium, per unit area of 
the plane boundary, is then 


j(T,»8) = ¥(T,.8) cos@, O<6< 7/2, (1) 


where 6 is the angle of incidence with respect to the 
normal to the plane. The current j has the same 
units as the fluence y. 


For the purposes of normalization, we define the 
integrated incident 41 fluence as* 


T 


4, = an [% if sinedd y(T,.6) 
0 0 


or 
One 4n Highs ¥(T,) R (2) 
fe) 


since y is assumed isotropic. 4, has the units cm °. 


Only particles with 0<6< 7/2 cross into a semi- 
infinite plane medium, so we define the integrated 
incident 2n fluence as 


1/2 ; 
on an f at, if sinede o(T,.8) 


a f at, ¥(T,) (3) 


and the total incident current as 


1/2 
J= an faty f sinede j(T,.8) 
) 0 
© 1/2 
= an f dT if sinédé cosé v(T_,8@) 
) re) 
fe) i) 
= aii dT, ¥(T,) (4) 
fe) 
Thus, 
0, = 25, = 4d (5) 


* 
A factor 2n comes from the integration over the 
azimuthal angle. 
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Results for Monoenergetic Radiation 


Electron Component 


Electron depth-dose curves and spectra were calcu- 
lated for incident energies Wye Wolls W525 W255 ila 4 
3, 5, 7, and 10 MeV, based on samples of 10,000 incident 
electrons in each case. The histories of all genera- 
tions of electrons were followed down to an energy of 
T,/32, at which point the deposition of the smal] 
amount of residual energy was treated approximately. 
rep histories were followed down to an energy of 

eV. 


In addition to the depth-dose distributions, the 
electron fluence spectra Fo(T,Z3T9), as a function of 
spectral energy T and depth z, set up within the 
medium were calculated for those electrons traveling in 
the forward direction (forward fluence denoted as Ffor) 
as well as for electrons traveling in all directions 
(total fluence denoted as Fiot). The data were obtained 
for depths z up to the electron mean range alg) 
the region in which the energy deposition is dominated 
by ionization (and excitation) losses of the primary 
electrons. Typical electron fluence spectra are shown 
in Fig. 1, for T, = 1 MeV, at both a depth near the 
entrance plane and a depth well into the medium. It can 
be seen that there is a sizeable backscattered component 
(the difference between the FEOt and the FLO curves) 
and a noticeable build-up of low-energy electrons. 

Using the notation of Jordan,’* we let Re(2. 1) 
be the aluminum dose as a function of depth z in the 
semi-infinite medium as directly calculated. The dose 
in thin tissue-equivalent detectors (very thin water 
layers) was obtained according to* ~~ 


H,0 
2 eae 
DE (2,T,) a D= (z,T,) 
T 
(0) H,0 
f reot(tyzsT,) L 2 (T.a) dT + N(z)A 
A 
x T > (6) 
0 
if Feotcr 7st) vA®(T,a) dT + N(z)A 
A e oO 


where L(T,A) is the electron restricted stopping power?? 
which includes only energy transfers less than A (in 
our case, A = 1,/32). In Eq. 6, N(z) is the number 

of electrons per unit depth falling below the cut-off 
energy A (the so-called track-ends). 


Similarly, the dose D_(z,To) 
finite slabs of thickness Zz 


at the edge of 
was approximated by 


** 
The mean range, also called the csda range because it 
is derived in the continuous-slowing-down approxi- 
mation, is defined as the integral of the reciprocal 
of the total stopping power over the energy range 
from the initial energy Tg down to some very small 
cut-off energy. Reference 11 contains an extensive 
tabulation of electron range values. 


kK 
This expression is closely related to the Spencer- 
Attix!? cavity-ionization theory. 


0 
f FLOT(T 237.) LAM(T,a)dT + N(2)A/2 
oM*(z,1,) = 0M (2,15) S-—--____—_ (7) 
@ 0 
If Feet (T25T.) LAL (T,a)dT + N(z)A 
A e {e} 
and by 
é 
) H,0 
Wie f eu epai) (h(i) clr & MrayVie 
LO Ota) (3) 


(Ore u 


: ALT A) dT + N(z)A 


In Eqs. 7 and 8, for the forward-moving electrons we use N(z)/2 in the track-end term since we find that at low 
spectral energies the fluence is nearly equally divided between the forward- and the backward-moving electrons. 


2/r, = 0-0.025 


z/r, = 0.425-0.45 


- ELECTRON FLUENCE F, /®,,, .MeV~' 


| 0.1 1 
| ELECTRON ENERGY T, MeV 


: 
FIG. 1. 


Electron spectra in a semi-infinite aluminum 
medium bombarded by a 1-MeV, isotropic 
electron beam. Calculated spectra Fe, as a 
function of electron spectral energy T, are 
shown for the depth intervals 0 < z < 0.025 ry 
and 0.425 <z<0.45 r, where r, is the 
electron mean range. At each depth, the lower 
curve is for electrons traveling in the 
forward direction, the upper curve includes 
electrons traveling in all directions. 


4898 


28) 


Bremsstrahlung Component, 


In order to obtain statistically significant results 
for the bremsstrahlung contribution to the dose, this 
contribution was not taken directly from the output of 
the ETRAN calculations, but indirectly in the following 
manner. A feature of the Monte Carlo code was exploited 
which allows for an arbitrary number of bremsstrahlung 
Photons per electron to be generated. This is compen- 
sated by the use of weight factors to correct back to the 
true production rate. In this way, separate calculations 
were made of the photon spectrum F.(k,23T9) as a 
function of photon spectral energy k and depth z for 
To = 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10 and 20 MeV, 
and for depths up to 30 g/cm?. These results are based 
on sample sizes of 10,000 incident electrons and from 
115,000 to 160,000 emitted bremsstrahlung photons. 


Figure 2 gives typical forward and total photon 
fluence spectra for T, = 1 MeV. At a shallow depth 
(~ 0.45 g/cm?, still within the electron range), soft 
bremsstrahlung photons are evident in the spectra with 
energies k < 10 keV. At larger depths, the spectrum 
hardens due to the strong photoelectric absorption of 
low-energy photons. 


The bremsstrahlung dose was then obtained from 
1) phe 

AX ut i tot 4 en 

D (z,T,) / Ee (k,25T,) k ae dk , 


co 


H50 


a 
H.0 ) Hee ik) 

2 sf tot ; wees 
LE Pali i ret (k,25T,) k EX dk, 


; 1, Ae 
QL for en 
D(z, : 023 , 
ise (ZiT) if B (DeBlg) ma Cl 
H0 
T 2 
H.0 ) Mae ES) 
2 z for : en 
p? (z,T,) il FF (k,25T,) k 2 — dk, (9) 


where en/p is the photon mass energy-absorption 
coefficient tabulated by Hubbell?*. 
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a 
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= 108 
F 1 10 100 1000 
@ INCIDENT PROTON ENERGY To, MeV 
3 . . . 
% 19 FIG. 3. Proton dose for fixed depths z in a semi- — 
FA infinite aluminum medium, plotted as a function 
Fs of incident energy To. The calculated results 
D. are given in terms of the dose per unit 2n 
isotropic incident fluence. 
10° 
~ 0.001 0.01 0.1 4 


PHOTON ENERGY k, MeV 


FIG. 2.  Bremsstrahlung spectra in a semi-infinite 
aluminum medium bombarded by a 1-MeV, isotropic 
electron beam. Calculated spectra Fy, asa 
function of bremsstrahlung photon spectral 
energy k, are shown for the depth intervals 
0.4<z<0.5 g/cm? and 9<z< 10 g/em?. At 
each depth, the lower curve is for photons 
traveling in the forward direction, the upper 
curve includes photons traveling in all 
directions. 


Basic Set of Depth-Dose Distributions 


Using the methods described above, a set of depth- 
dose distributions for aluminum shields were generated 
for monoenergetic sources. The distributions were 
smoothed in z and To. The overall uncertainty of the 
final results, including the uncertainties in cross 
sections used and those due to statistical fluctuations 
and smoothing of the Monte Carlo data, is estimated to 
be 10%. 


ELECTRON DOSE Da/®,,. rads (AL) cm2 


The aluminum dose expressed in rads*, as a function 
of To for fixed depths, is shown for a semi-infinite 
medium in Fig. 3 for protons and in Fig. 4 for elec- 
trons. The electron curves include contributions from 
primary and secondary knock-on electrons and from 
secondary electrons set in motion by bremsstrahlung 
photons. At z= 0, the proton dose reduces to the 
stopping power. However, for the semi-infinite medium, 0.1 1 10 


the electron dose at z = 0 is larger than the stopping 
power due to backscattering from greater depths. INCIDENT ELECTRON ENERGY. To Mev 


FIG. 4. Electron dose for fixed depths z ina semi- 
infinite aluminum medium, plotted as a function 


*n rad = 10° Gy = 1.602 x 10° MeV/g. The detector of incident energy To. The calculated results 
material will be specified within parentheses in the D= are given in terms of the dose per unit 2m 
unit of dose, e.g., rads (Ag). isotropic incident fluence. The dashed curve 

ae is the zero-depth dose based only on the 
In the straight-ahead approximation, D- = D- for stopping power and neglects backscattering 
protons. fe from greater depths. 


24 4899 


In Fig. 5, the electron dose as a function of To, 
at z= 0, 0.5 and 5 g/cm*, is compared for the cases 
of the aluminum dose in a semi-infinite aluminum medium, 
the aluminum dose at the edge of finite aluminum slabs, 
and the tissue dose in a semi-infinite aluminum medium. 
For very thin slabs the dose is less than the stopping 
power because of the escape of energy in the form of 
energetic knock-on electrons. We find that except for 
the region where the bremsstrahlung dominates, the 
tissue dose is roughly a factor 1.3 larger than the 
aluminum dose and that the dose for finite slabs is 
roughly 0.7 of that for a semi-infinite medium (there 
are variations with z and To of as muchas ~ + 20% 
from this approximate value for the D_/D, ratio). 
However, for the bremsstrahlung component, the situation 
is more complicated because of the strong variation with 
photon energy and detector material of the mass energy- 
absorption coefficient. 


107 


10° 


10° 
—Ds/®,,, . rads (AZ) 

: —--D—/®,,,. rads (AZ) 
oO 
4 ---0s/9,, . 
© 4910 rads (HO) 
Ww 
77) 
{e) 
fa 
5 11 
2 10 
= 
oO 
Ww 
a 
Ww 

10°12 

10°13 

10;ut 

INCIDENT ELECTRON ENERGY To. MeV 
FIG. 5. Comparison of electron dose distributions in 


semi-infinite aluminum media and at the edge 
of finite aluminum slabs, for both aluminum 
and tissue-equivalent detectors. Solid curves 
pertain to the aluminum dose D— at depths z 
in a semi-infinite aluminum medium; dot-dashed 
curves pertain to the aluminum dose D_ at 
the edge of aluminum slabs of thickness 2; 
dashed curves pertain to the tissue dose D. 
at depths z ina semi-infinite aluminum 
medium. The points for z= 0 are based only 
on the stopping power. 


* 
This factor can be estimated simply from stopping pewer 
tables. Based on the calculated fluence spectra, we 
find values from about 1.24 to 1.34. 
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Measurements were made by Van Gunten’® of the 
aluminum dose in hollow aluminum cubes, one face of 
which was irradiated by monoenergetic, isotropic electron 
beams. Those measurements made directly behind and at 
the center of the irradiated face would seem to closely 
correspond to D_, except possibly for a small additional 
contribution backscattered from the other sides.of the 
cube. The calculated results are compared in Table I to 
the properly normalized’® data of Van Gunten. For Ty = 1, 
2 and 3 MeV, the agreement is good, within the stated 
uncertainty of 10% for the experimental data. However, 
there is a discrepancy at 4 MeV where the measured doses 
are consistently lower than the calculated ones. We can 
find no ready explanation for this discrepancy. 


TABLE I. Comparison of calculated and measured dose 
behind aluminum slab shields irradiated by 
monoenergetic, isotropic electron beams. The 
experimental results of Van Gunten!® and the 
calculated results are given in units of 
10°” rads (AZ), and are normalized to a unit 
2m isotropic incident fluence. Results are 
given for various incident electron energies 
Ty and shield thickness z. 

To (MeV) z(mils*) Exptl. Dose Calc. Dose 

] 30 0.110 + 0.011 0.115 
] 60 0.00781 + 0.0078 0.00701 
2 30 0.164 + 0.016 0.167 
2 60 0.118 + 0.012 0.120 
2 125 0.0163 + 0.0016 0.0153 
3 30 0.169 + 0.017 0.181 
3 60 0.155 + 0.016 0.153 
3 125 0.0812 + 0.008) 0.0859 
4 30 0.149 + 0.015 0.189 
4 60 0.136 + 0.014 0.168 
4 125 0.0915 + 0.0092 0.122 
4 250 0.0146 + 0.0015 0.0252 

a : 2 ; 

1 mil = 0.00686 g/cm” of aluminum 


Calculations for Arbitrary Incident Spectra 


The monoenergetic-source results have been incorpor- 
ated into a computer code SHIELDOSE** which rapidly 
performs the necessary interpolation and integration for 
any incident spectrum. 


Slab Shields 


The results from the approximate proton calculations 
are compared in Fig. 6 with those calculated by Burrell?’ 
and by Alsmiller?® for aluminum slabs. The SHIELDOSE 
results are no more than 10-20% higher than these other 
results which include the effects of nuclear interactions. 
Such a discrepancy is consistent with the findings of 
Santoro et al’? who tested the effects of nuclear 
collisions and secondary particles in proton space- 
shielding calculations. 


Electron results for simple exponential incident 
spectra (To) = (1/2ma) exp(- To/a) are given in Fig. 7 
and are compared to those of Watts and Burrel]l*°. When 
the bremsstrahlung component is omitted the curves are 
generally in good agreement as should be expected, since 
both calculations are based on data obtained (indepen- 
dently) using ETRAN. In order to indicate the effect of 
bremsstrahlung on the depth-dose distribution, curves 
including this component are given also in Fig. 7. 


kk 
This code will be available from the Radiation Shielding 
Information Center, Oak Ridge National Laboratory, and 
from the National Space Science Data Center, Goddard 
Space Flight Center. 


PROTON DOSE Dg, rads (HO) 


ELECTRON DOSE 0z/,, , rads (H,O) cm2 


2. g/cm2 


Comparison of proton depth-dose distributions 
D. in aluminum slab shields. The dashed curve 
is from the present calculations, the solid 
curve from Alsmiller!®, and the points from 
Burrell’’. All results were calculated for an 
incident, jsotropic proton fluence (in 
EmisMeversnay) Of v(ilg)) = 2a5eX On ligt mak 
for 5 < Ty < 60 MeV and ¢(To) = 5.486 x 10°* 
peter for 60 < To < 1000 MeV, and pertain to 
irradiation from one side only. 


FIG. 


Previous bremsstrahlung calculations in space- 


shielding 
generally 
mation. 


problems have been done in various ways, but 
have been based on the straight-ahead approxi- 


As an example, we compare in Fig. 8 our results 


with those of Watts and Burrell*° for exponential 


incident electron spectra. 


number of 


distributions of the emitted bremsstrahlung and give 
various results ranging from those from their most 


elaborate 


estimate (based in part on our older Monte Carlo data” 
pertaining to the energy content of bremsstrahlung beams 
emerging from aluminum targets). 


These authors considered a 
models describing the energy and angular 


treatment to a simple order-of-magnitude 


As seen in Fig. 8, 


these results disagree with the SHIELDOSE results by as 


much as a 


their more approximate results fortuitously agree to 
within ~ 30%. 


Spherical Shields 


factor of two; however, it turns out that 


In application to spacecraft there is usually 


interest in geometries other than those of simple slab 


targets. 


beyond the scope of this work. 
been obtained for the dose at the center of uniform 
solid spheres using slab depth-dose distributions. 


Calculations for complicated geometries are 
However, results have 


This 


is accomplished using the relation 


The Electron dose D. as a function of depth 2z 
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a=4 MeV 


10° 


40°10 


0 1 2 3 4 5 6 7 
DEPTH z, g/cm2 


in 
a semi-infinite aluminum medium. Results are 
for exponential incident spectra, ¥(T) = const. 
x exp (- To/a), and are given in terms of the 
dose per unit 27 isotropic incident fluence. 
Dashed curves are from calculations of Watts 
and Burrell2° based on ETRAN depth-dose data; 
solid curves are from present calculations 
(upper curves include bremsstrahlung component, 
lower curves do not). 


an isotropic fluence of particles in the straight- -ahead 
approximation. 
merely takes into account that a sphere is irradiated 
from all directions while the semi-infinite medium is 
irradiated only from one side. 


The factor of 2 = $4,/¢2, in Eq. 10 


Proton doses at centers of spheres are shown in 
Fig. 9, where our results are compared with those 
calculated by Burrell?” and Hill? 
by ~ 10-20% for large radius 2z 
as in the slab case -- attributable to our neglect of 
attenuation due to nuclear interactions. 


The overestimate 
of our results, is -- 


In the case of electrons and bremsstrahlung, the 
use of Eq. 
more approximate. 
calculations specifically for the slab and spherical 
geometries in order to validate Eq. 
depths where the bremsstrahlung contribution is still 


10 is, because of multiple scattering effects, 
Jordan’? has made Monte Carlo 


10, at least up to 


' * 
Dee Dee) i d log D-(z') \ Eq. 10 is given by Jordan‘? in terms of the dose D_ 
° oo d log 2' ae A between two finite slab shields each of thickness 2, 
rather than in terms of the dose D— ina semi-infinite 
where D.(z) is the dose at the center of a sphere of medium. For protons in the straight- -ahead approximation, 
radius Zz. This relationship, which has been suggested one O53 Penge C ENON except for veny esmalliliig 
by Jordan!2 and references cited by him, is exact for where some differences in the backscattering becomes 
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evident, D_= 
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FIG. &. Bremsstrahlung dose D— as a function of 
depth z in * semi-infinite aluminum medium. 
Results are for exponential incident spectra 
with an e-folding energy a, and are given 
in terms of the dose per unit 2n isotropic 
incident fluence. Solid curves are from 
present calculations; broken curves are from 
the straight-ahead calculations of Watts and 
Burrel1?° (dashed curves from their most 
elaborate treatment, dot-dashed curves are 
their most approximate estimate). 


10 are his Monte Carlo results 
both for D- and D_ for the fission a 
v(To) = const x exp {-0.575 To -0.055 T Shown also 
are our results* based on our Monte Pune: calculations 
for D. and the use of Eq. 10 for D,. The disagreeme 
(~ 10%) in D, for small z is evidently due to our 
use of i rather than D_, while the disagreement 
(again ~ 0%) at large z for both D- and D 
apparently is the result of using different Monte Carlo 
calculations. 


small. Shown in Fig. 


The ratio of the sphere dose to the slab dose is 
plotted in Fig. 11 for the fission spectrum and for an 
exponential spectrum with a = 0.5 MeV. Shown also is 
the ratio derived from Jordan's data for the fission 
spectrum. We make the following observations pertinent 
to space-shielding calculations: (a) Using Eq. 10 for 
electrons appears to be an approximation good to within 
~ 10-20%. This is not to say that the straight-ahead 


Our dose in silicon was extimated from the dose in 
aluminum by multiplying by 1.03, which takes into 
account the electron stopping-power ratio. For the 
bremsstrahlung component also included in the results, 
this factor should actually be somewhat larger but has 
little effect on the comparison. 
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FIG. 9. Proton dose at the centers of aluminum spheres 


D,, as a function of sphere radius z. Dashed 
curves are from present calculations; solid 
curves are from those of Burrell!’; and 

points are from Hill??. Upper set of results 
are for an incident, isotropic proton fluence 
(in cm *MeV''sr'*) of (Tg) = 5.486 x 107% 

Ty °*°°3 the lower set are for v(Tg) = 1.257 x 
1022 To *°1%. Irradiation of the Sphere is 
assumed to be from all directions. 


approximation per se, i.e., neglecting angular deflec- 
tions, is accurate for describing electron penetration, 
but rather that, for the incident electron distributions 
of interest, electron penetration can be described by 
some functions (details of which we need not know) as 
if the electrons move along straight lines. (b) The 
argument in favor of using the results of Eq. 10 for 

the bremsstrahlung tail is the success of some form of 
the straight-ahead approximation in describing at least 
the shape of the depth-dose distribution as demonstrated, 
for example, in Fig. 8. (c) Whereas D ,/D5 is less 
than about 6 to 8 for protons, this ratio Ts as large 

as 13 for electrons and about 3 in the bremsstrahlung 
tail. 


Realistic Spectra 
Finally, we give an example of SHIELDOSE results 


for incident spectra kindly provided by E.G. Stassinopoulos 


of NASA/GSFC using the radiation-environment models and 
orbit-integration codes** developed at the National 
Space Science Data Center. Figure 12 shows the one- 
year aluminum dose, both D- and D,, for the omni- 
directional, orbit- - integrated, trapped electron fluence 
for a geostationary orbit with an altitude of 35790 km, 
an inclination of 0°, and for a parking longitude of 
160°W. The spectrum, shown in the insert of Fig. 12, 


FIG. 


ELECTRON DOSE , rads(Si) cm2 


10° 


SOLID SPHERE, 
D./ Ben 


SEMI-INFINITE SLAB, \ 
Os/ny 


Comparison of the electron dose in a semi- 
infinite medium (at depth z) and at the 
centers of spheres (radius z). Results are 
for aluminum targets (100 mils = 0.686 g/cm?) 
irradiated by an isotropic electron fluence 
having a fission spectrum, ¥(T,) = 0.057 
exp(-0.575 Ty - 0.055 Tp”) cm *MeV ‘sr *. 
Da'shed curves are Jordan's!* results from 
Monte Carlo calculations for each geometry. 
Solid curves are the present results based on 
Monte Carlo calculations for the semi-infinite 
medium and on the use of Eq. 10 for the 
spheres. 


was calculated for the epoch 1979.0 using the AEI7-HI 
model. 
period is (To) = 9.25 x 10° exp(- Ty/26.5 MeV) 

cm *MeV ‘sr ?, and corresponds to one anomalously 
large event predicted at the 95% confidence level. 


The 


The solar proton fluence for this orbit and 
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FISSION SPECTRUM 


100 200 300 400 500 600 700 800 $00 1000 


2, mils 


11. Ratio of the electron dose D, at centers of 
aluminum spheres (radius z) to twice the 
dose Di. in a semi-infinite aluminum medium 
(at depth z). The results, pertaining to 
aluminum detectors, are plotted out to large 
z where the bremsstrahlung dominates. Solid 
curves are our results based on Eq. 10 and 
are given both for incident electrons with a 
fission spectrum and for electrons with an 
exponential spectrum characterized by an 
e-folding energy a= 0.5 MeV. Dashed curve 
is the fission-spectrum results of Jordan}? 
who considered explicitly the two geometries. 
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INCIDENT ELECTRON FLUENCE giT,). cm=? MeV-' or~' 


10° 


O17 1 
T,. Mev 


BREMSSTRAHLUNG 
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12. Depth-dose distributions in aluminum targets 
for the radiation encountered during one year 
in a geostationary orbit. The orbit parameters 
and incident radiation spectra, are described 
in the text; the incident electron fluence is 
given in the insert. Results are given both 
for the dose D- at depth z ina semi- 
infinite medium and for the dose D, at the 
center of a sphere of radius z. The solar- 
proton dose, for 1 anomalously-large event, 
is given by the dashed curves. 
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Appendix B. Monoenergetic depth-dose data. 

Scaled, monoenergetic depth-dose distributions are given in the following 
tables. The results are normalized to a unit incident current (1 cm72) of 
particles crossing into the plane aluminum medium, and pertain to particles 
with an isotropic angular distribution incident from one side only. The dose, 
absorbed in thin detectors of various materials of interest as given both as 
a function of depth z in a semi-infinite medium, ptet(z), and as a 
function of thickness z at the transmission SuRRaCeiOr a finite thickness 
slab, Mea) Reference can be made to figure 1 for a further indication 
of irradiation and target geometry. The scaling is done in terms of the 
mean (c.s.d.a.) range els) in g/cm? (given in tables 1 and 2 of the 
main text) and incident energy, T. in MeV. The tables are organized as 


0) 
follows: 


Ag ‘ Ag 
Table B.1. eT) Bs (z/r,) = (e/a) De (z/r,) for protons. 
Units: dimensionless. 


H50 


H,0 
2 o 2 
Table B.2. (ey) Be (z/r,) = ra) D_ (z/r,) for protons. 


Units: dimensionless. 


Si a Si 
Table B.3. rs) D (z/r) = (r/T)) D~ (Zire) for protons. 


Units: dimensionless. 
$i0 
2/1) = Ca.) Di 


Units: dimensionless. 


Si, 
Table B.4. (r,/T)) Dine 


co 


(z/r,) for protons. 


Table B.5. Me /aley De (z/r,) for electrons. Units: dimensionless. 


co 


Table B.6. (ne) D*(z/r.,) for electrons. Units: dimensionless. 
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Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


Table 


B. 


B. 


B. 


B. 


B. 


B. 


B. 


B. 


Be 


B. 


10. 


Wc 


Wo 


20. 


: Oe le) a 


HA0 
(ro/T,) D_ 


(z/r,) for electrons. Units: dimensionless. 
50 
Weg) D_ (Zins) for electrons. Units: dimensionless. 
(nea) p>" (z/r ) for electrons. Units: dimensionless. 
(0) 4S(0) Fs (@) 
(r,/T,) D°"(z/r,) for electrons. Units: dimensionless. 
(2/1) for electrons. Units: dimensionless. 


Sid 
a ey De (z/r,) for electrons. Units: dimensionless. 


(UT) DY (z/T,) for bremstrahlung. Units: cm@/g. 


‘ i) Die (ey) for bremstrahlung. Units: cm/g. 


H50 


; rar) Die 27) for bremstrahlung. Units: cm@/g. 


oe) 


(1/T,) ACA) for bremstrahlung. Units: cm2/g. 


‘ (1/T,) p>"(z/T,) for bremstrahlung. Units: cm2/g. 


foe) 


; Ws) Deutz) for bremstrahlung. Units: cm@/g. 


Si 


0 
; (1/T,) D_ 2 (2/T,) for bremstrahlung. Units: cm@/g. 


(1/T,) D_ (z/T,) for bremstrahlung. Units: cm‘/g. 
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Appendix C. Sample Run 


]. Input Data 
Below is a listing of the input data for the run producing the results 
shown in figure 12 of Appendix A. Note that the monoenergetic depth-dose data 


set is read from unit 9. 


1 9 2 
42 2 

0.03000 0.04000 0.05000 0.06000 0.08000 0.10000 

0.15000 0.20000 0.30000 0.40000 0.50000 0.60000 

0.80000 1.00000 1.25000 1.50000 1.75000 2.00000 

2.50000 3.00000 3.50000 4.00000 4.50000 5.00000 

5.50000 6.00000 6.50000 7.00000 7.50000 8.00000 

8.50000 9.00000 9.50000 10.00000 12.50000 15.00000 

17.50000 20.00000 22.50000 2? yeas 27.50000 30.00000 

2.000 1000.000 2.000 1000.000 30 0.100 7.000 101 

GEOSTAT, Sa 88 KM, INCL=0,PLONG=160W; SP: 1AL(95%) » TP: NONE, EL: AEI7-HIC79) 

3 28 1.00000E 03 3.15360E 07 4 
0.0000E 00 0.0000E 00 0.0000E 00 | 
2.6500E 01 2.4500E 10 0.0000E 00 
1.0000E-01 2.0000E-01 3.0000E-01 4.0000E-01 5.0000E-01 6.0000E-01 i 
7.0000E-01 8.0000E-01 9.0000E-01 1.0000E 00 1.2500E 00 1.5000E 00 
1.7500E 00 2.0000E 00 2.2500E 00 2.5000E 00 2.7500E 00 3.0000E 00 
3.2500E 00 3.5000E 00 3.7500E 00 4.0000E 00 4.2500E 00 4.5000E 00 
4.7500E 00 5.0000E 00 5.5000E 00 6.0000E 00 6.5000E 00 7.0000E 00 
7.5343E 04 5.7226E 04 3.5768E 04 2.0637E 04 1.20846E 04 7.7592E 03 
5.3738E 03 4.0677E 03 3.3050E 03 2.0840E 03 9.4707E 02 4.4939E 02 
2.2097E 02 1.0415E 02 5.1592E 01 2.6350E 01 %1.2667E 01 5.7739E 00 
2.9428E 00 1.8846E 00 1.4374E 00 1.1547E 00 8.4948E-01 6.2276E-01 
4.8761E-01 3.4227E-01 2.1946E-01 1.2189E-02 0.0 0.0 


2. Output Data 
The output listing, the significant portions of which are given below, is 
mostly self-explanatory. In addition to monitoring the input data, the total 
flux and the average energy is given, both for the protons and electron 


incident fluxes: 
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Sample output listing: 


EAV = f 
- 


min 


IMAX IUNT 
42 2 
SHIELD DEPTH (G/CM2) 
0.03000 0.04000 
0.15000 0.20000 
0.80000 1.00000 
2.50000 3.00000 
5.50000 6.00000 
8.50000 9.00000 
17.50000 20.00000 
| EMINS EMAXS EMINP 
2.000 1000.000 2.000 
2.000 DIO. 9)9) 2.000 


+ 


max 
INT SPEC = if ley) Clits 
T 


-05000 
30000 
-25000 
.50000 
-50000 
.50000 
-50000 


NWAw-oo 


EMAXP 
1000.000 
999.999 


GEOSTAT,35790 KM, INCL=0,PLONG=160W; 


JSMAX JPMAX JEM 
3 0 


E(MEV) 
0.0 0.0 0.0 
SOLAR PROTON SPECTRUM (/ENERGY/CM2) 
| 2.6500E#01 2.4500E+10 0.0 
INT SPEC — EAV(MEV) 
2.2719E+10  28.50006 
E(MEV) 
1.0000E-01 2.0000E-01 3.0000E-01 
7.0000E-01 8.0000E-01 9.0000E-01 
1.7500E+00 2.0000E+00 2.2500E+00 
3.2500E+00 3.5000E+00 3.7500E+00 
4.7500E+00 5.0000E+00 5.5000E+00 
| ELECTRON SPECTRUM (/ENERGY/CM2/TIME) 
_—-7.5343E+04 5.7226E+04 3.5768E+04 
5.3738E+03 4.0677E+03 3.3050E+03 
2.2097E+02 1.0413E+02 5.1592E+01 
2.9428E+00 1.8846E+00 1.4374E+00 
4.8761E-01 3.4227E-01 2.1946E-01 
INT SPEC —_ EAV( MEV) 
1.9162E+07 0.34290 


AX 


EUNIT 


min 


0.06000 0.08000 
0.40000 0.50000 
1.50000 1.75000 
4.00000 4.50000 
7.00000 7.50000 
10.00000 12.50000 
25.00000 27.50000 
NPTSP EMINE 
301 0.100 7.000 
301 0.100 7.000 


28 1.00000E+03 3.15360E+07 


PnM— LL 


——NWNNP 


-0000E-01 
-9000E+00 
-5000E+00 
-O000E+00 
-0000E+00 


-0637E+04 
-0840Et03 
-6350E+01 
- 1547E+00 
-2189E-02 
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ANU 


Co o— 


-0000E-01 
-2500E+00 
-7500E+00 
-2500E+00 


-2084E+04 
-4707E+02 
-2667E+01 
-4948E-01 


RUKH 


AufL~ 


Tinax 
iL Van (To) qT, 
Thin 


0.10000 
0.60000 
2.00000 
5.00000 
8.00000 
5.00000 
0.00000 


EMAXE NPTSE 


101 
101 


SP: 1AL(095%),TP: NONE, EL: AEI7-HIC79) 
TINTER 


-0000E-01 
-5000E+00 
-0000E+00 
-5000E+00 


-7592E+03 
-F939E+02 
-7739Et00 
-2276E-01 


ADJUSTED VALUES 


OO+O29E°S 
00+H9E8°9 
O0+HLG4°8 
LO+adbel ol 
LO+adescb 
LO+d106°1L 
LO+d0£S 2 
LO+429b°¢ 
LO+4S26°% 
LO#d005°S 
LO+d0L8°S 
LO+adb15°9 
LO+tGZ1L0°2 
LOtadtliod 
LO+td25S 8 
LO+tH01S 6 
20+4990°1 
20+d602 1 
COtALLE 1 
20+498S 1 
20+A8S8° 1 
cOtdole'?¢ 
c0+d289 2 
CO+AG8E £ 
COtHLG9° 4 
20+H692°S 
20+d29G 2 
E£O+AGEL Ll 
£0+d262°2 
£0+E8S5°S 
5b0+H629° 1 
DO+HD6L ES 
bOtHLC> 9 
GO+dGld 1 
GO+46f>°¢ 
GOtGZ12L°S 
90+4590°1 
90+4G55° 1 
90+4890°¢ 
90+406GS°2 
90+4261°¢ 
90tH9E1°4 


IVLOL 


00+4560°S 
00+4055°9 
OO+FHZEL 8 
LO+agco'l 
Lo+d00e | 
LO+4859° 1 
LO+dzol 2 
10+4602°2 
LO+d0GS'°¢ 
LO+td8GZ°¢ 
LO+d£86°E 
lo+dlee >» 
LOtdf65° 4 
LO+d 582° 
LO+adb0L°S 
10+46495°S 
LO+4dSS8°S 
LO+4dc0£°9 
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Appendix D. Listing of the computer code. 


SHIELDOSE, 12 FEB 80. 
DEVEEORED BY S. Mi SELTZER: 


IDET = 1, AL DETECTOR 

2, H2-O DETECTOR 

3, SI DETECTOR 

4, SI-O02 DETECTOR 
INCIDENT OMNIDIRECTIONAL FLUX IN ZENERGY/CM2/UNIT TIME 
CSQLAR-FLARE FLUX IN /ENERGY/CM2). 
EUNIT IS CONVERSION FACTOR FROM /ENERGY TO /MEV, 

Ex Ge seUNIn e— NO00 SE IREUX TS KEV. 

TINTER IS MISSION DURATION IN MULTIPLES OF UNIT TIME. 


NATIONAL BUREAU OF STANDARDS. 


DIMENSION EP(C30),RPC30),DP(30,51) 


ZP(51),ER(80),REC80),EEC10), 


DE(10,41,2),2E(41),2B(60),EB(10),DB(10,60,2),2M(50),2(50), 
TPLC301),TPC301),RINPC301),TELC101),TEC101),RINEC101),S0L(301), 
EPS€101),5¢€101),SPG(301),60301),SEG(101),D0S0L(50,2), 
DOSP(50,2),DIN(60,301),GP(301,50),GE(101,50,2),G6B(101,50,2), 
DOSE(50,2,2),D0SB(50,2,2),ZMMN(50),Z2L(50),GARB(60) 


CALL ERRSET (208,256,-1,1) 

ZCON=0.001*2.540005*2.70 

ZMCON=10.072.70 

RADCON=1.6021892E-08 

PRINT 10 

FORMAT (€°0 IDET NUT IPRNT') 

READ 30, IDET,NUT,IPRNT 

PRINT 30, IDET,NUT,IPRNT 

IDET1L=IDET+1 

REWIND NUT 

PRINT 20 

FORMAT C'OMPMAX LPMAX KMAX MEMAX LEMAX MBMAX LBMAX"*) 
READ CNUT,30) MPMAX,LPMAX,KMAX,MEMAX, LEMAX,MBMAX, LBMAX 
PRINT 30, MPMAX,LPMAX,KMAX,MEMAX, LEMAX,MBMAX, LBMAX 
FORMAT (1216) 

PRINT 40 

FORMAT C'OPROTON DOSE DATA CENERGIES INCREASING)'") 
PRINT 50 

FORMAT (€'0 ECMEV) 

DO 110 M=1,MPMAX 

READ (NUT,60) EPCM),RPCM) 
PRINT 60, EPC(M),RPCM) 
FORMAT C1P6E12.4) 
EP(M)=ALOGCEP(M) ) 
RPCM)=ALOGC(RP(M) ) 

DO 70 ID=1,IDET 

READ (NUT,60) CDP(M,L),L=1,LPMAX) 


RCG/CM2)") 


CONTINUE 

GO TO (80,90), IPRNT 

PRINT 60, (DPC(M,L),L=1,LPMAX) 
UG GLEE Gha¢) GO TO 20 


DO 100 ID=IDET1,4% 

READ (NUT,60) CGARBC(L),L=1,LPMAX) 
CONTINUE 

CONTINUE 
ZPFAC=1.0/FLOATCLPMAX-1) 

DO 120 L=1,LPMAX 
ZP(L)=FLOATCL-1)*ZPFAC 

PRINT 130 

FORMAT C'OQELECTRON DOSE DATA’) 
PRINT 140 

FORMAT C'OECMEV)') 

READ (NUT,150) CERCK),K=1,KMAX) 
PRINT 150, CERCK),K=1,KMAX) 
FORMAT (6F12.5) 

PRINT 160 

FORMAT C'ORANGECG/CM2)") 

READ (NUT,60) CRECK),K=1,KMAX) 
PRINT 60, CRECK),K=1,KMAX) 
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170 
180 


190 
200 
210 
220 


230 
240 


250 
260 
270 
280 
290 


300 


310 
320 
330 


340 
350 


360 


370 
380 


390 


400 
410 
420 


430 


440 
450 


DO 170 K=1,KMAX 

ERCK) =ALOGCER(K)) 

RECK) =ALOGCRECK) ) 

PRINT 180 

FORMAT C'OECMEV), THEN ROXD/E CENERGIES INCREASING) ') 
DO 270 M=1,MEMAX 

READ (NUT,150) EECM) 

PRINT 150, EECM) 

EECM) =ALOGCEE(M) ) 

DO 210 ID=1,IDET 

DO 200 N=1,2 
READ (NUT,190) C(CDECM,L,N),L=1,LEMAX) 
FORMAT C1P8E10.3) 
CONTINUE 
CONTINUE 
GO TO (220,240), IPRNT 

DO 230 N=1,2 

FRINT 190, (DECM,L,N),L=1,LEMAX) 

CONTINUE 

LF CLUDE GES) CO TOe270 

DO 260 ID=IDET1,4 

DO 250 N=1,2 

READ CNUT,190) CGARB(L),L=1,LEMAX) 

CONTINUE 

CONTINUE 

CONTINUE 

ZEFAC=1.0/FLOATCLEMAX-1) 

DO 280 L=1,LEMAX 

ZECLJ=FLOATCL-1)*ZEFAC 

PRINT 290 

FORMAT C*OBREMSSTRAHLUNG DOSE DATA CENERGIES INCREASING) ') 
PRINT 300 

FORMAT C'O02Z2/7E (G/CM2 KEV)'") 

READ (NUT,310) (ZB(L),L=1,LBMAX) 

PRINT 310, (ZBCL),L=1,LBMAX) 

FORMAT C1P8E9.2) 

DO 320 L=1,LBMAX 

ZBCL)=ALOG(1000.0X%ZBC(L)) 

PRINT 330 

FORMAT C'OECMEV), THEN D/E (CM27G)") 

DO 420 M=1,MBMAX 

READ CNUT,150) EBCM),RENORM 

PRINT 150, EBCM),RENORM 

EBCM)=ALOGCEBC(M) ) 

DO 350 ID=1,IDET 

DO 340 N=1,2 

READ (NUT,310) (CDB(M,L,N},L=1,LBMAX) 

CONTINUE 

CONTINUE 

GO TO ¢€360,380), IPRNT 

DO 370 N=1,2 

PRINT 310, (DBCM,L,N),L=1,LBMAX) 

CONTINUE 

DO 390 N=1,2 

DO 390 L=1,LBMAX 

DB(M,L,N)=ALOGCRENORM¥DBC(M,L,N)) 

LE CLDER.Gi.4)) GO) 10.420 

DO 410 ID=IDET1,4 

DO 400 N=1,2 

READ CNUT,310) CGARBC(L),L=1,LBMAX) 

CONTINUE 

CONTINUE 

CONTINUE 

PRINT 430 

FORMAT ('0 IMAX IUNT") 

READ 30, IMAX,IUNT 

PRINT 30, IMAX,IUNT 

GO TO (440,470,500), IUNT 

PRINT 450 

FORMAT C"OSHIELD DEPTH (MILS)") 
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460 


470 
480 


490 


500 
510 


520 
530 
540 


550 


570 


580 


590 


600 
610 
620 


READ 150, (ZMCI), 

PRINT 150, (ZMCI) 

DO 460 I=1,IMAX 

Z(I)=ZCONXZMC( I) 

ZMNCI)=ZC1)*ZMCON 

GO TO 530 

PRINT 480 

FORMAT C'OSHIELD DEPTH (G/CM2)') 

READ 150, (2¢1I),I=1,IMAX) 

PRINT 150, (ZC(1I),I=1,IMAX) 

DO 490 I=1,IMAX 

ZMCI)=Z¢01I)7ZCON 

ZMMCI)=Z01)*ZMCON 

GO TO 530 

PRINT 510 

FORMAT C'QSHIELD DEPTH (MM)"') 

READ 150, (ZMMCI),I=1,IMAX) 

PRINT 150, (ZMMCI),I=1,IMAX) 

DO 520 I=1,IMAX 

ZCI =ZMMCID/ZMCON 

2ZMNC1I)=Z201)7ZCON 

DQ 540 I=l1,IMAX 

ZLCID=ALOG(ZC1)) 

PRINT 550 

FORMAT ¢'O EMINS EMAXS EMINP EMAXP NPTSP EMINE 
ENAXE NETSE*) 

READ 560, EMINS,EMAXS,EMINP, EMAXP,NPTSP,EMINE,EMAXE,NPTSE 

PRINT 560, EMINS,EMAXS,EMINP, EMAXP,NPTSP,EMINE, EMAXE,NPTSE 

FORMAT (4F10.3,16,2F10.3,16) 

EMINU=AMINICEMINP, EMINS ) 

ENAXU=AMAX1CEMAXP, EMAXS ) 

DEP=ALOGCEMAXU/EMINU)/FLOATCNPTSP-1) 

NFSTS=ALOGCEMINS/ZEMINU)/DEP+0.5 

NFSTS=NFSTS#1 

NLSTS=ALOGCEMAXS/EMINU)/DEP+0.5 

NLSTS=NLSTS#1 

NLENS=NLSTS-NFSTS+1 

HFSTP=ALOGCEMINP/EMINU)/DEPt0.5 

NFSTP=NFSTP#+1 

NLSTP=ALOGC EMAXP/EMINU)/DEP+0.5 

NLSTP=NLSTP+1 

NLENP=NLSTP-NFSTPtHI1 

EMINUL=ALOGCEMINU) 

DELP=DEP/73.0 

ICALL=1 

DO 570 NP=1,NPTSP 

TPLONP)=EMINUL+FLOATCNP-1)*DEP 

TPOCNP)=EXPCTPLONP)) 

CALL SPOL CTPLCNP),EP,RP,MPMAX,ICALL,ANS) 

RINPCNP)=EXPCANS) 

PRINT 580, TPCNFSTS),TPCNLSTS),TPCNFSTP),TPCNLSTP),NPTSP,EMINE, EMA 


T=1,IMAX) 
»T=1,IMAX) 


IXE,NPTSE 


FORMAT (4F10.3,16,2F10.3,16,"' ADJUSTED VALUES") 
DO 620 L=1,LPMAX 

ICALL=1 

DO 610 NP=1,NPTSP 

IF (TPLONP).LT.EPCMPMAX)) GO TO 590 
DINCL,NP)=DPCMPMAX,L) 

GO TO 610 

IF CTPLONP).GT.EP(1)) GO TO 600 
DINCL,NP)=DP(1,L) 

GO TO 610 

CALL SPOL CTPLCNP),EP,DP(1,L),MPMAX,ICALL,DINCL,NP)) 
CONTINUE 

CONTINUE 

DO 660 NP=1,NPTSP 

ICALL=1 

DO 650 I=1,IMAX 

ZRIN=ZCI)/RINPCNP) 

IF (ZRIN.LT.1.0) GO TO 640 
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1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 

1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 


630 
640 


650 
660 


680 


SMI 
eHow 
ooo 


720 


730 


740 
750 


760 


770 
780 
790 


800 
810 
820 
830 
840 


850 


860 


1 
PRINT 870 


GPCNP,1I)=0.0 

GO TO 650 

CALL SPOL (ZRIN,ZP,DIN(C1,NP),LPMAX,ICALL,ANS) 
IF CANS.LT.0.0) GO TO 630 
GPCNP,I)=TPCNP)XANS/RINPCNP) 

CONTINUE 

CONTINUE 

EMINEL=ALOGCEMINE) 
DEE=CALOGCEMAXE)-EMINEL)/FLOATCNPTSE-1) 
DELE=DEE/3.0 

ICALL=1 

DO 670 NE=1,NPTSE 

TELCNE) =EMINEL+FLOATCNE-1)*DEE 

TECNE) =EXPCTELCNE) ) 

CALL SPOL CTELCNE),ER,RE,KMAX,ICALL,ANS) 
RINECNE) =EXP CANS) 

DO 820 N=1,2 

DO 710 L=1,LEMAX 

ICALL=1 

DO 700 NE=1,NPTSE 

IF (TELCNE).LT.EECMEMAX)) GO TO 680 
DINCL,NE)=DECMEMAX,L,N) 

GO TO 700 

TE CREUCNED Gi EEC!) G0 1:0) 690 
DINCL,NE)=DEC1,L,N) 

GO TO 700 

CALL SPOL (TELCNE),EE,DEC1,L,N),MEMAX,ICALL,DINC(L,NE)D) 
CONTINUE 

CONTINUE 

DO 750 NE=1,NPTSE 

ICALL=1 

DO 740 I=1,IMAX 

ZRIN=ZCI)/RINECNE) 

TE CZRING E05 GO 0) 27:30 
GECNE,1I,N)=0.0 

GO TO 740 

CALL SPOL (ZRIN,ZE,DINC1,NE),LEMAX, ICALL,ANS) 
TRSGANS ETO 0D GOO 720 
GECNE,I,N)=TECNE) XANS/RINECNE) 

CONTINUE 

CONTINUE 

DO 790 L=1,LBMAX 

ICALL=1 

DO 780 NE=1,NPTSE < 

IF (TELCNE).LT.EBCMBMAX)) GO TO 760 
DINCL,NE)=DBCMBMAX,L,N) 

GO TO 780 

TE ChrEECNEX SGT SEBCID) GO 10-7770 
DINCL,NE)=DBC1,L,N) 

GO TO 780 

CALL SPOL (TELC(NE),EB,DBC(1,L,N),MBMAX,ICALL,DIN(L,NE)) 
CONTINUE 

CONTINUE 

DO 810 NE=1,NPTSE 

ICALL=1 

DO 800 I=1,IMAX 

ZBIN=ALOG(ZCI)/TECNE) ) 

CALL SPOL (ZBIN,ZB,DIN(1,NE),LBMAX,ICALL,ANS) 
GBCNE,I,N)=TECNE)XEXPCANS) 

CONTINUE 

CONTINUE 

PRINT 840 

FORMAT Ctr?) 

PRINT 850 

FORMAT (€'0") 

READ (5,860, END=1420) 

PRINT 860 

FORMAT (72H 
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870 FORMAT C'OJSMAX JPMAX JEMAX EUNIT TINTER’) 
READ 880, JSMAX,JPMAX, JEMAX, EUNIT, TINTER 
PRINT 880, JSMAX,JPMAX, JEMAX, EUNIT,TINTER 
880 FORMAT (316,1P2E12.5) 
IF (TINTER.LE.0.0) TINTER=1.0 
DELTAS=RADCONXDELP/74.0 
DELTAP=TINTERXRADCONXDELP/4.0 
DELTAE=TINTERXRADCONXDELE/4.0 
IF CEUNIT.LE.0.0) EUNIT=1.0 
ISOL=2 
IF (JSMAX.LT.3) GO TO 900 
ISOL=1 
PRINT 140 
READ 60, CEPS(J),J=1,JSMAX) 
PRINT 60, CEPS(J ypeleatky JSMAX) 
PRINT 890 
890 FORMAT C'OSOLAR PROTON SPECTRUM (/ENERGY/CM2)'") 
READ 60, (S(J),J=1,JSMAX) 
PRINT 60, (SC€J),J=1,JSMAX) 
CALL SPECTR (JSMAX,EPS,S,EUNIT,DELP,NLENS,TPCNFS7S),TPLONFSTS),SOL 
LONFSTS)) 
900 ITRP=2 
IF (JPMAX.LT.3) GO TO 920 
ITRP=1 
PRINT 140 
READ 60, CEPS(J),J=1,JPMAX) 
PRINT 60, CEPSCJ),J=1,JPMAX) 
PRINT 910 
910 FORMAT C'OTRAPPED PROTON SPECTRUM (/ENERGY/CM2/TIME)") 
READ 60, (S(€J),J=1,JPMAX) 
PRINT 60, (S(J),J=1,JPMAX) 
CALL SPECTR (CJPMAX,EPS,S,EUNIT,DELP,NLENP,TPCNFSTP),TPLONFSTP),SPG 
1CNFSTP)) 
920 ILEC=2 
IF (JEMAX.LT.3) GO TO 940 
ILEC=1 
PRINT 140 
READ 60, CE 
PRINT 60, (¢ 
PRINT 930 
930 FORMAT ('0 
READ 60, ¢ 
PRINT 60, 
CALL SPECT 
940 GO TO (980, 
950 DO 960 NP=N 
960 SOLCNP)=0.0 
DO 970 J=l, 
1,1IMA Ae 
)= 
0 


970 DOSOL(CI,J 
GO TO 101 
980 DO 1000 I=1,IMAX 
DO 990 NP=NFSTS,NLSTS 
990 GCNP)=SOLCNP)XGPCNP,I) 
CALL INT (DELTAS, GCNFSTS), NLENS,DOSOL(I,1)) 
1000 CONTINUE 
CALL SPHERE (2L, DOSOL(1, 1), IMAX, DOSOL(1,.2)) 
10190 GO TO (1050, 1020), ITRP 
1020 DO 1030 NP=NFSTP,NLSTP 
1030 SPGCNP)=0.0 
DO 1040 J=1,2 
DO 1040 I=1,IMAX 
1040 DOSP(I,J)=0.0 
GO TO 1080 
1050 DO 1070 I=1,IMAX 
DO 1060 NP=NFSTP,NLSTP 
1060 GCNP)=SPGCNP)XGPCNP,I) 
CALL INT CDELTAP,GCNFSTP),NLENP,DOSP(I,1)) 
1070 CONTINUE 
CALL SPHERE (ZL,DOSP(1,1),IMAX,DOSP(1,2)) 


» JEMAX) 
1, JEMAX) 


Li ed 


E 
S 
( 
R 
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2780 
2790 
2800 
2810 
2820 
23830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3230 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 


1080 GO TO (1110,1090), ILEC A 3480 
1090 DO 1100 J=1,2 A 3490 
DO 1100 N=1,2 A 3500 

DO 1100 I=1,IMAX A 3519 
DOSECI,N,J)=0.0 A 3520 

1100 DOSBCI,N,J)=0.0 A 3530 
GO TO 1160 A 3540 

1110 DO 1150 N=1,2 A 3550 
DO 1130 I=1,IMAX A 3560 

DO 1120 NE=1,NPTSE A 3570 
GCNE)=SEGCNE)*GECNE,I,N) A 3580 

1120 SPGCNE)=SEGC(NE)*XGBCNE,I,N) A 3590 
CALL INT (DELTAE,G,NPTSE,DOSECI,N,1)) A 3600 

CALL INT (DELTAE,SPG,NPTSE,DOSBC(I,N,1)) A 3610 

1130 CONTINUE A 3620 
GO TO (1150,1140), N A 3630 

1140 CALL SPHERE (ZL,DOSEC(1,N,1),IMAX,DOSE(1,N,2)) A 3640 
CALL SPHERE (2ZL,DOSB(1,N,1),IMAX,DOSB(1,N,2)) A 3650 

1150 CONTINUE A 3660 
1160 J=1 A 3670 
DO 1340 N=1,2 A 3689 

GO TO (1170,1190), N A 3690 

1170 PRINT 1180 A 3700 
1180 FORMAT C'1DOSE AT TRANSMISSION SURFACE OF FINITE ALUMINUM SLAB SHI A 3710 
LELDS") A 3720 

60 TO 1210 A 3730 

1190 PRINT 1200 A 3740 
1200 FORMAT C'1DOSE IN SEMI-INFINITE ALUMINUM MEDIUM") A 3750 
1210 GO TO (1220,1240,1260,1280), IDET A 3760 
1220 PRINT 1230 A 3770 
1230 FORMAT ¢C'ORADS AL®™) A 3780 
GO TO 1300 A 3790 

1240 PRINT 1250 A 3800 
1250 FORMAT C*ORADS H20"') A 3810 
GO TO 1300 A 3820 

1260 FRINT 1270 A 3830 
1270 FORMAT C'ORADS SI") A 3840 
GO TO 1300 A 3850 

1280 PRINT 1290 A 3860 
1290 FORMAT C'ORADS SI-02') A 3870 
1300 PRINT 1310 A 3880 
1310 FORMAT ('0 ZC(MILS) Z2¢(MM) Z(G/CM2) ELECTRON BREMS A 3890 
12 J EESBR TRP PROT SOL PROT EL+BR+TRP TOTAL") A 3900 
PRINT 850 A 3910 

DO 1330 I=1,IMAX A 3920 
DOSEB=DOSECI,N,J)+DOSBCI,N,J) A 3930 
DOSEBP=DOSEB+DOSP(I,J) A 3940 
DOST=DOSEBP+DOSOL(I,J) A 3950 
PRINT 1320, ZMCI),Z2MMC(I),Z(1),DOSECI,N,J),DOSBCI,N,J),DOSEB,DOSP(I A 3960 
1,J),DOSOL(I,J),DOSEBP,DOST A 3970 
1320 FORMAT COP3F11.3,1P7E11.3) A 3980 
1330 CONTINUE A 3990 
1340 CONTINUE A 4000 
J=2 A 4010 

N=2 A 4020 
PRINT 1350 A 4030 

1350 FORMAT ('1172 DOSE AT CENTER OF ALUMINUM SPHERES") A 4040 
GO TO (1360,1370,1380,1390), IDET A 4050 

1360 PRINT 1230 A 4060 
GO TO 1400 A 4070 

1370 PRINT 1250 A 4080 
GO TO 1400 A 4090 

1380 PRINT 1270 A 4100 
GO TO 1400 A 4110 

1390 PRINT 1290 A 4120 
1400 FRINT 1310 A 4130 
PRINT 850 A 4140 

DO 1410 I=1,IMAX A 4150 
DOSEB=DOSECI,N,J)+DOSBCI,N,J) A 4160 
DOSEBP=DOSEB+DOSP(I,J) A 4170 
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1410 
1420 


10 

20 

30 
40 


50 


: 60 
70 


} 10 


20 


_ 


DOST=DOSEBP+DOSOL(I,J) 


PRINT 1320, 2MC(I),ZMM(I),Z(1),D0SECI,N,J),DOSBCI,N,J),DOSEB,DOSPC(I 


»J),DOSOLC(I,J), DOSEBP, DOST 
CONTINUE 

GO TO 830 

STOP 

END 


SUBROUTINE SPECTR (JMAX,EPS,S,EUNIT,DELTA,NPTS,T,TL,SP), 
SUBROUTINE SPECTR CJMAX,EPS,S,EUNIT,DELTA,NPTS,T,1TL,SP) 


DIMENSION: EPS(1),S¢€1),T¢1),TLO1),SP(1),60301) 
TESGERSCT)).Gi-0.0) GO TO 20 
ALPHA=S(1) 

BETA=S(2) 

IF (BETA.LE.0.0) BETA=1.0 
BETA=BETA/ALPHA 

DO 10 N=1,NPTS 

SPCN)J=TCN) XBETAXEXP(-TCN)/ALPHA) 
GCN) =TCN)¥SPCN) 

GO TO 50 

DO 30 J=1,JMAX 

EPSCJ)=ALOGCEPS(J)) 
SCJI=ALOGCEUNITXS(J)) 

ICALL=1 

DO 40 N=1,NPTS 

CALL SPOL C(CTLC(N),EPS,S,JMAX,ICALL,ANS) 
SPOON) =TOCN)*EXPCANS) 

GCND=TCN)IXSPCN) 

CALL INT CDELTA,SP,NPTS,SIN) 

CALL INT (DELTA,G,NPTS, EBAR) 
EBAR=EBAR/SIN 

PRINT 60 

FORMAT ('0 INT SPEC EAVCMEV)") 
PRINT 70, SIN,EBAR 

FORMAT C1PE12.4,0PF12.5) 

RETURN 

END 


SUBROUTINE SPOL (S,X,Y,N,IN,T), 15 JAN 71. 
SUBROUTINE SPOL (S,X,Y,N,IN,T) 


CUBIC SPLINE INTERPOLATION WITH PARABOLIC RUNOUT. 


DIMENSION X(1),YC1),E°€101),UC101) 
GO TO (10,50), IN 
IN=2 


B2=X(Jt1)- 
C2=(Y¥CJt1)- 
B=X(J+1)-X¢ 
D=(C2-C1)/B 
C=B1/B 
B1=B2 

C1=C2 
P=CXECJ-1)+2.0 
ECJ)=(€C-1.0)/P 
UCJI=CD-CXUCJ-1))7P 
ECN) =UCN1)7(1.0-ECN1)) 
DO 30 KK=1,N1 

K=N-KK 


30 ECK)=ECK) XECK+1)+UCK) 
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IF (X€1).GT.X€N)) GO TO 40 


IDIR=0 280 
MLB=9 290 
MUB=N 300 
GO TO 50 310 
40 IDIR=1 320 
MLB=N 330 
MUB=0 340 


50 IF (S.GE.X(MUB+IDIR)) GO TO 390 


IF (S.LE.XC(MLBt+1-IDIR)) GO TO 100 360 
ML=MLB 370 
MU=MUB 

GO TO 70 390 


60 IF CIABSCMU-ML).LE.1) GO TO 110 
70 MAV=CMLTtMU)Z2 
IF (S.LT.XCMAV)) GO TO 80 


ML=MAV 430 
GO TO 60 440 
80 MU=MAV 450 
GO TO 60 460 
90 MU=MUB+2XIDIR 470 
GO TO 120 480 


ololkekelolelelololeleleloleleleleolelelelelene) 
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100 MU=MLBt+2xC1-IDIR) 


SUBROUTINE INTCDELTA,G,N,RESULT), (CN=1). 
SUBROUTINE INT (DELTA,G,N,RESULT) 


DIMENSION G(1) 30 
NLI=N-1 40 
NL2=N-2 50 


IF CFLOATCN)-2.0XFLOATCN/2)) 90,90,10 
10 IF (N-1) 20,20,30 
20 SIGMA=0.0 

GO TO 80 
30 IF (N-3) 40,40,50 


40 SIGMA=G(1)+4.0*G(2)+G6(3) 110 
GO TO 80 120 
50 SUM4¢=0.0 130 


DO 60 K=2,NL1,2 

60 SUM4=SUM4+G(K) 
SUM2=0.0 
DO 70 K=3,NL2,2 

70 SUM2=SUM2+G(K) 
SIGMA=G(1)4+4.0*XSUM4¢4+2.0*XSUM2+G(N) 

80 RESULT=DELTAXSIGMA 
RETURN 

90 IF CN-2) 100,100,110 

100 SIGMA=1.5*(G(1)+G6(2)) 
GO TO 80 

110 IF (N-4) 120,120,130 

120 SIGMA=1.125*(GC1)+3.0%G(2)+3.0*%G(3)+6(4)) 
GO TO 80 

130 IF (N-6) 140,140,150 

140 EU Slate gaa colette ache ea aa edad Cala abot AG 
y O 80 

150 IF (N-8) 160,160,170 

160 SIGMA=G(1)+3.875*G(2)42.625*G6(3)4+2.625%6(4)43.875*%G(5)42.0%G(6)44. 


10*GC7)+G6(8) 330 
GO TO 80 340 
170 SIG6=G(1)+3.875*%G(2)+2.625%G(3)+2.625%6(4)+3.875*%G(5)+G6(6) 350 
SUN4=0.0 360 


DO 180 K=7,NL1,2 
180 SUM4=SUM4+G(K) 

SUM2=0.0 

DO 190 K=8,NL2,2 
190 SUM2=SUM2+G(K) 


wieleiveiwleivelviolwelivivioivlolvivoiolveliviclololojololeloleleleloleleoleloelelololelojojel~) 
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SIGMA=SIG6+G6(6)+4.0XSUMN4+2.0*XSUM2+G(N) 420 
GO TO 80 430 
END 440 
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SUBROUTINE SPHERE (ZL,DOSE,IMAX,DOSPH), 30 AUG 79 
SUBROUTINE SPHERE (ZL,DOSE,IMAX,DOSPH) 

DIMENSION 2ZL(1),D0SE(1),D0SL(50),DERV(50),DOSPH(1) 
DO 10 I=1,IMAX 

HEP CDOSE CT) SEES 0.0) (GO T0)-20 
DOSLCI)=ALOGCDOSECI)) 

T=IMAXt+1 

IMIX=I-1 

te GUikx is) GO 10) 40 

CALL SPLDRV (ZL,DOSL,DERV,IMIX) 

DO 30 I=1,IMIX 

DOSPHCI)=DOSECI)*(1.0-DERVC(I)) 

IMIX1L=IMIX+1 

IF CIMIX1.GT.IMAX) RETURN 

DO 50 I=IMIX1,IMAX 

DOSPHC(T)=0.0 

RETURN 

END 


SUBROUTINE SPLDRV (X,Y,U,N), 28 AUG 79. 

SUBROUTINE SPLDRV (X,Y,U,N) 

CUBIC SPLINE WITH PARABOLIC RUNOUT. 
U CONTAINS DERIVATIVES AT KNOTS 

DIMENSION X(1),YC1),E(101),UC1) 

NI=N-1 

EC1)=1.0 

UC1)=0.0 

BIl=X(2)-X(1) 

C1l=CY(C2)-YC1))7B1 

DO 10 J=2,N1 

B2=X(Jt1)-XCJ) 

C2=CYCJ+t1)-Y¥CII)7B2 

B=XCJt1)-XCJ-1) 

D=(C2-C1)/B 

C=B1/B 

B1=B2 

C1=C2 

P=CXECJ-1)+2.0 

ECJ)=C€C-1.0)7P 

UCJ)=CD-CXUCJ-1))7P 

ECN)=UCN1)/7(01.0-ECN1)) 

DO 20 KK=1,N1 

K=N-KK 

ECK)=ECK) XECK+1)4UCK) 

B2=X(K+1)-X(K) 

UCKI=ECYCK+1LI-YCK) )7B2-B2*%(2. OXECKIFECKH1)) 

CONTINUE 

UCNI=CYCNI-YONT) 7B2+B2*%02.0XECN)J+ECNI) D 

RETURN 

END 
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